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Classical Hamiltonian systems with a mixed phase space and some asymmetry may exhibit 
chaotic ratchet effects. The most significant such effect is a directed momentum current or 
acceleration. In known model systems, this effect may arise only for sufficiently strong chaos. In 
this paper, a Hamiltonian ratchet accelerator is introduced, featuring a momentum current for 
arbitrarily weak chaos. The system is a realistic generalized kicked rotor and is exactly solvable 
to some extent, leading to analytical expressions for the momentum current. While this current 
arises also for relatively strong chaos, the maximal current is shown to occur, at least in one case, 
precisely in a limit of arbitrarily weak chaos. 



PACS numbers: 05.45.-a, 05.45.Mt, 05.60.-k, 45.05. 



I. INTRODUCTION 

Classical [H-Q and quantum 0, ^-[13 Hamiltonian 
ratchets have attracted a considerable theoretical interest 
during the last decade. Also, several kinds of quantum 
ratchets have been experimentally realized using atom- 
optics methods with cold atoms or Bose-Einstein con- 
densates [l3i - ll6j |. The classical Hamiltonian ratchet ef- 
fect is a directed current in the chaotic region generated 
by an unbiased force (having zero mean in space and/or 
time) and due to some spatial and/or temporal asym- 
metry [lref. This is analogous to the ordinary ratchet 
effect [13, [3 > but with deterministic chaos replacing the 
usual noisy environment. Dissipation, an important in- 
gredient in ordinary ratchets for breaking time-inversion 
symmetry, is absent in Hamiltonian ratchets. 

A well studied class of systems are those described by 
time-periodic Hamiltonians H(x,p,t) for which both the 
force F = —dH/dx and the velocity v = dH/dp are 
periodic in x and F has zero mean over (x,t). The 
classical ratchet current is usually defined as the aver- 
age of v over (x,p,t), where (xjA is restricted to the 
chaotic region, see, e.g., Refs. 0,9- ^ is assumed that v 
is bounded, e.g., by Kolmogorov-Arnol'd-Moser (KAM) 
tori. Then, necessary conditions for the ratchet current 
to be nonzero are the breaking of some symmetry and 
a mixed phase space featuring "transporting" stability 
islands which propagate in the x direction [3j. In the 
presence of bounding KAM tori, one can get, in princi- 
ple, a well-defined ratchet current also in near-integrable 
regimes, corresponding to relatively weak and local chaos. 

A different and much more significant Hamiltonian 
ratchet effect was discovered in work [6] for generalized 
kicked-rotor systems satisfying the well-known KAM sce- 
nario. Namely, for sufficiently strong kicking, there exist 
no KAM tori bounding the chaotic motion in the momen- 
tum (p) direction, and one then gets strong global chaos. 
In addition, there may arise transporting "accelerator- 
mode" islands [19] propagating in the p direction. This 
can lead, under some asymmetry conditions, to a "ratchet 
acceleration", i.e., a nonzero mean momentum velocity 



(rather than the usual position velocity v) of the global 
chaotic region @ (see more details in Sec. II). Quantum 
analogs of the classical ratchet acceleration were found 
in several systems either for special, "quantum- 

resonance" values of a scaled Planck constant H @-[Hi or 
for generic values of h [12{ . Quantum-resonance ratchet 
accelerators have been experimentally realized in recent 
works [H.[l5l|. 

In this paper, we show for the first time that the phe- 
nomenon of ratchet acceleration is not limited to strong- 
chaos regimes. We introduce a realistic Hamiltonian 
system exhibiting this phenomenon most significantly in 
near-integrable regimes, corresponding now to arbitrar- 
ily weak but global chaos. The system is a generalized 
kicked rotor whose force function has zero mean and is 
characterized by two nonintegrability parameters b\ and 
£>2 • A global chaotic region in the p direction arises also 
for arbitrarily small values of these parameters, i.e., the 
KAM scenario is not satisfied. As b\, 62 — > 0, this 
"non-KAM" [2(| system tends to the well-known ellip- 
tic sawtooth map [111, [2^, which has been used as a 
paradigmatic model of "pseudochaos" (dynamical com- 
plexity with zero Lyapunov exponent) [22, 23] in studies 
of both classical [22|, [24j and quantum [25| systems. We 
show that accelerator-mode islands exist for arbitrarily 
small b\ and 62- Then, when the system is asymmet- 
ric (b\ 62), a ratchet acceleration A may arise for ar- 
bitrarily weak chaos. In one particular case, we derive 
analytical expressions for A as a function of b\ and 62- 
Paths of maximal A in the (b\,b-i) parameter space are 
determined. We then show that, in sharp contrast with 
the systems considered in work p|, A is most significant 
for relatively small Lyapunov exponent and that its max- 
imal value is attained precisely in a limit b\, 62 - > of 
arbitrarily weak chaos. 

This paper is organized as follows. In Sec. II, we 
give a short background on Hamiltonian ratchet accel- 
erators. In Sec. Ill, we introduce our general model 
system and describe its basic properties. In particular, 
in Sec. IIIC we derive the existence conditions for the 
main accelerator-mode islands of the system. In Sec. IV, 
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analytical expressions for the ratchet acceleration A in 
one case are obtained for all values of the parameters. 
In Sec. V, we show that the maximal value of A is at- 
tained in a limit of arbitrarily weak chaos. Conclusions 
are presented in Sec. VI. Detailed derivations of several 
analytical results are given in the Appendices. 

II. BACKGROUND ON HAMILTONIAN 
RATCHET ACCELERATORS 

The concept of Hamiltonian ratchet accelerator was 
introduced in Ref. [6|by adaptation of a formalism de- 
veloped in Refs. [2, 3]. We give here a self-contained 
summary of these works, leading to the main result ([6]) 
below, a sum rule for the ratchet acceleration A. We 
shall focus on realistic models, the generalized kicked- 
rotor systems with scaled Hamiltonian H = p 2 /2 + 
KV{x, t) Yl^L-oo — s )' wnere K is the nonintegra- 
bility parameter and the potential V(x,t) is periodic in 
x, V(x + l,t) = V(x,t). Particular cases of these sys- 
tems where considered in Ref. [6]. The map for H from 
t = s — Otoi = s+1 — Ois given by 

M : p s+1 = p s +Kf s {x s ), x s+ i = x s +p s+ i mod(l), (1) 

where the force function f s (x) = —dV(x,t = s)/dx. Be- 
cause of the periodicity of V(x, t) in x, f s (x) satisfies the 

l 

ratchet (zero- flux) condition: (f s (x)) — J f s (x)dx = 0. 

As in Ref. @ , we shall assume that the kicking parameter 
K is large enough that all the rotational ("horizontal") 
KAM tori are broken. Then, there are no barriers to mo- 
tion in the p direction, leading to a global and strongly 
chaotic region. These barriers cannot exist if there are 
"accelerator modes", i.e., orbits that are periodic under 
map ^ in the following sense: 



Ps+m — Ps 



Xs+r, 



(2) 



where m is the period and w, the winding number, is 
an nonzero integer. Periodic orbits can be defined in the 
generalized way © (w / 0) due to the obvious period- 
icity of the map ([1} in p with period 1. If an accelerator 
mode is linearly stable, each of its points {x s ,p s ) will be 
usually surrounded by an island T s , an "accelerator- mode 
island" (AI). Because of ©, l s + m is just I s translated 
by w in the p direction. For arbitrary initial conditions 
zo = (xq,po) in phase space, the mean acceleration (mo- 
mentum current /velocity) in n iterations of ([TJ is 



A„(z ) = 



Pn - Pa 



n 



(3) 



and the average of ([3]) in some region 1Z with area S-r. is 

1 ' ■ ■ ■■ (4, 



(A n )n = — / A n (z )dz , 
&n J 
n 



In the case that TZ is an AI I with winding number w, it 
follows from Eqs. ©-0) that 



T , A \ W 

lim {A n ) i = v = — . 

n— >oo m 



Now, because of the periodicity of the map ([TJ in p with 
period 1, one can take also p s +i modulo 1 in ([1}, leading 
to a map M on the unit torus T 2 : < x, p < 1; this is 
the unit cell of periodicity of map ([1]). The reduced phase 
space T 2 can be fully partitioned into the global chaotic 
region C with area Sc and all the stability islands 1^ 
with areas Sj, where j labels the island: Sc + Y] a Sj = 1. 

The average acceleration ([5]) ofl^' is i/j = Wj/rrij, where 
Vj = for a normal (non-accelerating) island. We then 
have the following sum rule relating Vj to the ratchet 
acceleration A — {A)c — lim n _ HX) (A n )c of the global 
chaotic region: 



(6) 



Eq. ([B]) is easily derived from the obvious relation 
{A n ) T 2 = S c {A n ) c + J2j S j ( A n)iU) by taking n -t oo 
and using {A n )j2 = 0, a result following straightfor- 
wardly from the map ([!]): 



\s=l / J2 s=0^ 2 

= - / d*.Ux.)=Q, (7) 

71 n J 



where we used area preservation (dzo = dz s ), the in- 
variance of T 2 under M, and the ratchet condition 
(fs{x)) = 0. An immediate consequence of relation ([6]) 
is that A vanishes if the map (|T|) is invariant under in- 
version, (x,p) — > {—x,—p), i.e., one has the inversion 
(anti) symmetry f s (—x) = —f s (x). This is because under 
this symmetry for each AI with mean acceleration Vj ^ 
there exists an AI with the same area but with mean ac- 



celeration 



As we shall see in the next sections for 



(5) 



a simple case of f s {x), A is generally nonzero when AIs 
are present and inversion symmetry is absent. 

III. THE GENERAL MODEL SYSTEM AND 
ITS BASIC PROPERTIES 

A. General 

The general model system introduced and studied in 
this paper is the generalized kicked-rotor system de- 
scribed by a simple map ([1]): 

M : p s+ i = p s +Kf(x s ), x s+ i = x s +p s+ i mod(l), (8) 

where < K < 4 and, for < x < 1, 

l\X for < x < &i, 

f(x) — ^ c - x for b± < x < 1 - b 2 , (9) 
l 2 {x - 1) for 1 — 6 2 < x < 1, 

with f(x + l) = f(x). Here b\ and b 2 are positive param- 
eters with bi+b 2 < 1 while l\, l 2 , and c, also positive, are 
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fixed by requiring f(x) to be continuous and to satisfy 
1 

the ratchet condition J f{x)dx =0 (see Appendix A): 

o 

(l-fri)(l-6i-fr 2 ) (1 - 6 2 )(1 - 6i - 6 2 ) 

bl{2 - bl -b 2 ) ' h b 2 {2-b x -b 2 ) ' 



Kicked systems with a smooth piecewise linear force 
function such as dU have been studied either on the 
phase plane [26[ or on a cylindrical phase space [27| . 
corresponding to the very special case of map (J8J) with 
b\ = b 2 = 1/4. Apparently, however, these systems have 
not been considered yet in the context of Hamiltonian 
ratchet transport, i.e., for general values of b\ and b 2 
with b\ ^ b 2l leading to an asymmetric force function 
©■ This general system is realistic since it may be ex- 
perimentally realized using, e.g., optical analogs as pro- 
posed in Ref. [26]. As we shall see below, the system 
generally does not satisfy the KAM scenario assumed in 
Sec. II, i.e., it is a non-KAM system. 

B. Phase space and limit cases 

The phase space of the map (JSj) in the basic periodicity 
torus T 2 (0 < x, p < 1) is illustrated in Fig. 1 for some 
values of the parameters. We clearly see in all cases a con- 
nected chaotic region encircling T 2 in both the x and p 
directions, implying global chaos and the nonexistence of 
KAM tori bounding p. An understanding of this numeri- 
cal observation will be achieved here and in Sec. IIIC. We 
first consider here the map ([5]) in the limit of b\, b 2 — > 0. 
From Eqs. (|10|) one has l\b\, l 2 b 2 , c — > 1/2 in this limit, 
so that the function © tends to the sawtooth 

f{x) = 1/2 -x (0<x<l), f(x + l)=.f(x), (11) 

with discontinuity at x — 0. The map © with (fTTj) 
and < K < 4 is the well-known elliptic sawtooth map 
(ESM) [H H [13, [13 having the property that its lin- 
earization DM is a constant 2x2 matrix with eigenvalues 
A± on the unit circle: 

A± = exp(±ia), 2 cos(a) = 2 - K. (12) 

This means that orbits of the ESM which do not cross 
the discontinuity line x — lie on ellipses with average 
rotation angle a. In general, however, an orbit will cross 
the x = line. Then, the combination of the mod(l) 
operation in © with the local ellipticity of the ESM will 
usually lead to a complex dynamics with zero Lyapunov 
exponent, known as "pseudochaos" (22[. The phase space 
generally consists of the pseudochaotic region, associated 
with all iterates of the discontinuity line [2lJ, and a set 
of islands. More specifically, one has to distinguish be- 
tween three main cases of the ESM, illustrated in Fig. 
2 for the same values of K as in Fig. 1: (a) The "in- 
tegrable" case of integer K = 1, 2, 3 [corresponding to 
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FIG. 1: Global chaotic regions of the map ([8]) within the unit 
torus of periodicity < x, p < 1 for bi = 0.05, b 2 = 0.02, 
and three values of K corresponding to the following values 
of a in Eq. G2]): (a) a = 2n/3 (K = 3); (b) a = 6tt/5 
(K « 3.618); (c) a = n(VE-l)/2 (K w 2.7247). These values 
of q represent the three main cases discussed in the text. The 
"left" (V) and "right" (R) period-1 accelerator-mode islands 
(AIs, see Sec. IIIC) are indicated in each case. 



a/(27r) = 1/6, 1/4, 1/3 in (|12j> ]: in this case, no pseu- 
dochaos arises and the phase space consists just of a finite 
number of "separatrix" lines (iterates of the discontinuity 
line) bounding a finite number of islands, see Fig. 2(a). 
(b) The case of non-integer K with rational a/(2n) in 
(IT21 ; in this case, numerical work [2l[ indicates that one 
has an infinite set of islands and that the pseudochaotic 
region is a fractal with zero area, see Fig. 2(b) and ex- 
act results for the fractal dimension of such regions in 
other maps with discontinuities [23] . (c) The case of ir- 
rational a/(27r); here one typically has again an infinite 
set of islands but the pseudochaotic region appears nu- 
merically to cover a finite area [2l[, see Fig. 2(c). Since 
the momentum p assumes all values on the discontinu- 
ity line and is thus unbounded, the pseudochaos [or the 
separatrix in case (a)] is global. 

For finite and small b\ and b 2 , the continuous map 
© may be considered as a perturbed ESM, with the 
discontinuity line replaced by a vertical strip B of width 
b\ + b 2 in T 2 (see also caption of Fig. 2) : 

B : <p < 1 and < x < h or l-b 2 < x < 1. (13) 
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FIG. 2: Global "pseudochaotic" regions of the ESM for the 
same values of a (or K) as in Fig. 1. In practice, these regions 
were generated by iterating a large initial ensemble using the 
map (JHJ with very small b\ —b^ — 10~ 7 . The initial ensemble 
covered uniformly the vertical "hyperbolic" strip (|13[) . The 
left and right AIs are again indicated. 



This should be contrasted with the perturbed ESM in 
Ref. [24| for which the discontinuity line is not removed 
by the perturbation. The linearization DM of ([5]) is again 
a constant 2x2 matrix in each of of the three intervals 
in Eq. (|9]). In the middle interval, it is the same matrix 
as for the ESM, with stability eigenvalues (fT2j) . In the 
other two intervals, where the strip (|T5]) is located, DM 
can be easily shown to have real positive eigenvalues A± 
with, say, A+ > 1 and A_ = A^ 1 < 1, i.e., there is local 
hypcrbolicity. One can then expect that already for small 
b\ and b 2 a global chaotic region with positive Lyapunov 
exponent will emerge from the vertical strip (|13[) and 
will replace the global pseudochaos (or separatrix) for 
bi = b 2 = 0. This can be clearly seen by comparing Figs. 
1 and 2. The nature of the chaotic region will be further 
discussed in the next sections, where it will be shown 
numerically that the Lyapunov exponent indeed tends to 
zero as b\, b 2 —> 0. 

C. Accelerator-mode islands (AIs) and their 
existence conditions 

We show here that AIs for the map {5} rigorously exist 
in broad ranges of the parameters, including arbitrarily 



small values of b\ and b 2 . This exactly implies global and 
arbitrarily weak chaos. We shall consider only period-1 
AIs, associated with stable accelerator modes satisfying 
Eq. © with m = 1 and w ^ 0. As we shall see, there 
appear to be no higher-period AIs at least in the case of 
K = 3 on which we shall focus from next section on. The 
initial conditions (iro,j?o) for m = 1 stable periodic orbits 
in Eq. @ must necessarily lie in the middle interval in 
Eq. {9J, &i < xo < 1 — b%, since only in this interval the 
matrix DM exhibits stability eigenvalues (fT2"|) . Then, 
from Eqs. ©, ©, and ©, we get: 



w 

x =c-—, p = 0mod(l), 
K 



h < c- — < 1 - b 2 . 
K 



(14) 



(15) 



For w — one has a non-accelerating stable fixed 
point (xo — c,po = 0), the center of a normal (non- 
accelerating) island. Let us show that w may take only 
two nonzero values and this only in some interval of K: 



w 



±1, 2<K<4. 



(16) 



In fact, from Eqs. (|9|) and (|T0|) it follows that the max- 
imal value of \f(x)\ is max(Zi&i, l 2 b 2 ) < 1/2. Then, 
since w — Kf(xo) from Eqs. ((2|) and ([8]), we have 
I to I < [-^72], where [ ] denotes integer part. This im- 
plies, for < K < 4, that w may take the only nonzero 
values of ±1 provided 2 < K < 4. 

Now, according to Eq. (1T41 for xq, the values of w = 1 
and w = — 1 should correspond, respectively, to a "left" 
(L) and "right" (R) AI, see Figs. 1 and 2. An explicit 
existence condition for the left AI (w — 1) is derived, 
after some simple algebra, from the left inequality in f| 15[) 
using Eq. (fTUl) for c: 



b 2 < F(h) EE M. 



{l-2K)b x +K -2 



K -I- Kbi 



(17) 



see also note [28|. It is easily verified that the right in- 
equality in (fT5j) is identically satisfied. Similarly, the ex- 
istence condition for the right AI is b\ < F(b 2 ). One thus 
has three cases (compare with Fig. 3 for K = 3): 

(a) Both AIs L and R exist (see, e.g., Figs. 1 and 2) if 



h < F(bi) and b x < F(b 2 ). 



(18) 



Clearly, this will be always satisfied for K > 2 and suffi- 
ciently small bi and b 2 since F(bi) w (K — 2)/(K— 1) for 
&i -C 1 in Eq. (fTT|); thus, both AIs exist in the arbitrarily 
weak chaos regime. For K = 3, this case corresponds to 
the domain LR in Fig. 3. 

(b) Only one AI, say the right one R, exists (as, e.g., 
in Fig. 4) if 

b 2 > F(h) and h < F(b 2 ). (19) 
Similarly, if b 2 < F(bi) and b x > F(b 2 ) only the left AI 
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reads 6i < F(b\), which can be significantly simplified: 
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FIG. 3: The curves b 2 = F(bi) and 6i = F(b 2 ) [with F(fci) 
given by Eq. (|17[) for AT = 3], defining the domains of exis- 
tence of the left and right period- 1 AI for K = 3 in the (6i, 62) 
plane. In domain L (R), only the left (right) AI exists. In 
domain LR, both AIs exist. No AIs exist elsewhere. 



1 1 

h< 2-K- 



(22) 



It follows from condition (|2"2"|) that no period- 1 AIs can 
exist if bi = b 2 > 1/4, for any value of K in the relevant 
interval of 2 < K < 4. This is consistent with the known 
fact that bounding KAM tori exist for some K if &i = 
62 = 1/4 (2?| , which is apparently the only case of map 
([8]) studied until now. 

IV. RATCHET ACCELERATION FOR K = 3 

In this section, the ratchet acceleration A in the case 
of K = 3 will be calculated analytically in the framework 
of a plausible assumption (see below) , supported by ex- 
tensive numerical evidence and exact results. To use the 
sum rule ([S]), we first identify the global chaotic region 
C in the basic periodicity torus T 2 . Let us denote by C 
the set of all iterates of the vertical strip B in Eq. iflF 
under M, i.e., the map (|SJ) modulo T 2 (see Sec. II): 




FIG. 4: Global chaotic region for K = 3, 61 = 0.226, and 
62 = 0.02. For these values of 61 and 62 only the right period- 
1 AI exists, compare with Fig. 3. 



L exists. For K = 3, this case corresponds to the domain 
R or L in Fig. 3. 
(c) No AIs exist if 



62 > F(bi) and h > F(b 2 ). 



(20) 



It is easy to show from the expression for F(pi) in Eq. 
(QUI (see also note [13) that F(6i) < (K-2)/(K-l)-bi- 
One then gets from Eqs. (fT8]) - (|20[) a simple necessary 
condition for the existence of at least one AI: 



h + b 2 < 



K -2 
K-l 



(21) 



For the symmetric system (b\ — b 2 ), the condition p7|) 



C 



u 



M S B. 



(23) 



Exact results for the set (j23[) are derived in Appendices 
B-E. Here we note that orbits which never visit B (and 
thus also C) are all stable since they lie entirely within 
the middle interval in ([9]) where the linearized map DM 
has stability eigenvalues (fL2|) . Thus, the global chaotic 
region C must be entirely contained within C, in agree- 
ment with our expectation at the end of Sec. IIIB. Our 
extensive numerical studies indicate that C is indistin- 
guishable from C, compare, e.g., Figs. 1(a) and 4 with 
Figs. 11 and 12 in Appendix C. In fact, finite-time Lya- 
punov exponents of orbits starting from initial conditions 
covering B uniformly were found to be all strictly pos- 
itive. We shall therefore assume in what follows that C 
precisely coincides with C. The rest of phase space out- 
side C consists of no more than three stability regions 
[see, e.g., Figs. 1(a) and 4]: The left AI L (w = 1), the 
right AI R (w = — 1), and a normal island (w = 0) lying 
between L and R. Using the sum rule © with Vj = Wj 
[since m — 1 in Eq. ([5])]. we then get a formula for the 
ratchet acceleration: 



.4 



St 



Sc 



(24) 



Exact expressions for the areas Sl, Sr, and Sc are de- 
rived in Appendices D and E using simple geometry; see 
Eqs. ijiS]). ([50]) . and (JUJ-flU) there. Inserting these ex- 
pressions in formula (|2~i]) . we obtain after some algebra 
explicit results for A in different cases: 
(a) If both AIs exist, i.e., case (TT8|) . 



.4 



(b 1 -b 2 )[l-3(6 1 +b 2 )] 
2(2-6 1 -6 2 )(6i+6 2 ) • 



(25) 
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(b) If only one AI, say the right one, exists, i.e., case 



.4 



(2 - 36 2 - 3c) 2 



6(6i 



6(&i 



b 2 y 



(3c -36! - l) 2 



(26) 



(c) If no AIs exist, i.e., case (|20l) . A — 0, of course. 

In general, the results (|25|) and (|26| were found to agree 
very well with numerical calculations of A, see examples 
at the end of the next section. This is additional evidence 
for the validity of the basic assumption above concerning 
the chaotic region, C = C. 

V. MAXIMAL RATCHET ACCELERATION 
FOR ARBITRARILY WEAK CHAOS 

In this section, we show that the maximal ratchet ac- 
celeration A for K — 3 is attained in a limit bx, b 2 — >• 
of arbitrarily weak chaos. In Fig. 5, we plot |A| as 
function of bx and 6 2 using formulas (f25j) . (|26|) . and 
A(b 2l bi) = —A{bx,b 2 )- The Lyapunov exponent a of 
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FIG. 5: (Color online) Pseudocolor plot of as function of 
b\ and 62 for K = 3. The thin solid lines (defining the three 
domains L, R, and LR) are the same as those in Fig. 3. The 
thick solid line in domain LR is the "maximal" path ()27|) on 
which A is given by Eq. (|28[l . The dashed line is the maximal 
path ^2(^1) [defined similarly to (|27[) ] on which A < 0. 

the chaotic region as function of 61 and b 2 was calculated 
numerically with high accuracy and is plotted in Fig. 6. 
As one could expect, a vanishes in the limit of 61, b 2 — > 0, 
where the map flS) tends to the ESM (see Sec. IIIB). It is 
clear from Fig. 5 that \ A\ assumes its largest values in the 
parameter domain LR, where both AIs exist. We shall 
therefore focus on this domain in which A is given by 
Eq. (j2"5")l . We shall hrst calculate analytically the value 



FIG. 6: (Color online) Pseudocolor plot of the Lyapunov ex- 
ponent a of the chaotic region as function of 61 and 62 for 
K = 3. 



of b\ where |^4(6i,&2)| is maximal at fixed b 2 \ this will 
define a path 61(62) in the (61,62) plane [a path 62(61) 
can be similarly defined]. We then show that |A(6i,&2)| 
is maximal on this path in the limit of 61, 62 — > 0. Let 
us take the partial derivative of the function (l2~5j) with 
respect to 61 and require that dA/dbi = 0. After a te- 
dious but straightforward calculation, we find that the 
latter equation reduces to a quadratic one with the only 
positive root: 



61(62) 



2 [56 2 (1 - 6 2 )] 1/2 -62(7-662) 
5-66, 



(27) 



The path (f27|) corresponds to the lower curve in Fig. 5, 
with 61 > 62. This curve starts at 61 = 62 = 0, with 61 w 
2^/62/5 for &2 1, and terminates at 61 = 6 2 = 1/6, 
on the boundary of the LR domain. For 62 < 1/6, we 
find that d 2 A/db\ < at the value ^ of 61, which thus 
corresponds to a local maximum. From Eqs. (|25[) and 
(|2"T)) . the ratchet acceleration on the path ([2"Tf is: 



A(b 2 ) 



2(J 



{5-6[56 2 (l-6 2 )] 1/2 } 2 
{5 - 462 - [56 2 (1 - 6 2 )] 1/2 } 



(28) 



In the limit of 62 — > (61 
(EH): 



2^62/5), we get from Eq. 



lim A{b 2 ) = 1/4. 



(29) 



After a simple but lengthy calculation, we find that the 
function ^ satisfies dA/db 2 < for 6 2 < 1/6. Thus, 
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A(p2) decreases monotonically from 1/4 (at b\ = b 2 = 0) 
to (at 61 = b 2 = 1/6) on the path (j2"T|) . Since this path 
gives the single extremum (a local maximum) of A(b±, b 2 ) 
for bi > 62 at fixed 62 and since A(bi, b 2 ) = for 61 = b 2 , 
we conclude that in the lower part (61 > 62) of the LR 
domain ^4(61,62) > and ^(61,62) assumes its maximal 
value of 1/4 in the limit 61, b 2 — > of arbitrarily weak 
chaos on the path (f!?7"l) . 

Since ^4(&2,6i) = —^(61,62), in the upper (b 2 > bi) 
part of the LR domain A(bi,b 2 ) < and A(bi,b 2 ) as- 
sumes its maximal negative value of —1/4 in the limit of 
61, b 2 — > on a path 62(61) (the upper curve in Fig. 5), 
defined similarly to 61(62)- The difference between the 
limiting values of A(bi,b 2 ) on the two paths reflects the 
discontinuity of the ESM, i.e., the map ([8]) in the limit 
of 61, 62 — > 0. In general, \A(bi,b 2 )\ can assume in this 
limit all values < 1/4 on other, "non- maximal" paths. 
For example, on the straight-line path 61 = 62/a, where 
a is some arbitrary constant, we find from Eq. (|25p that 



i- An \ (1 - a ) 

km A(b 2 ) = —, - 

b2 _H) v ; 4(1 



a) 



(30) 



We remark that the path (j2"T)) is tangent to the 61 axis at 
b\ = b 2 = 0, since 61 « 2y/b 2 /5 for b 2 <C 1. Similarly, the 
second maximal path 62(61) is tangent to the 62 axis in 
this limit. Thus, as expected, the maximal value of \ A\ — 
1 /4 is associated with the largest possible "asymmetry" , 
61/62 = 00 or 62/61 = 00 [a — or a — 00 in Eq. (|30)]. 

Figs. 7 and 8 show plots of A versus the Lyapunov 
exponent a for small 62 on both the maximal path (|2T[) 
and the path 61 = 36 2 . We see in both plots an excellent 
agreement between the values of A calculated numerically 
and those calculated from formulas (|2"5j) . (|26[) . and 



VI. CONCLUSIONS 

In this paper, we have introduced a realistic non-KAM 
system exhibiting, in weak-chaos regimes, the most sig- 
nificant Hamiltonian ratchet effect of directed accelera- 
tion. The system, defined by the generalized standard 
map © with ©, may be viewed as a perturbed ellip- 
tic sawtooth map (ESM) with a perturbation that re- 
moves the ESM discontinuity. Then, the global weak 
chaos featured by the system may be generally consid- 
ered as a perturbed global pseudochaos. Our main exact 
result is that for K = 3 the maximal ratchet acceler- 
ation A is attained precisely in a limit b\, b 2 — > of 
arbitrarily weak chaos with infinite asymmetry parame- 
ter (61/62 = 00 or 62/61 = 00). Despite this fact, the 
limiting system is interestingly the completely symmet- 
ric ESM (see phase spaces in Fig. 2). By continuity 
considerations, one expects that at least for values of K 
sufficiently close to K — 3 one should again observe a 
significant increase of the absolute value \A\ of the accel- 
eration as the chaos strength decreases. We have verified 
this numerically in parameter regimes where good ac- 
curacy could be achieved within the limitations of our 
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FIG. 7: Circles: Numerical results for A versus the Lyapunov 
exponent a on the maximal path (|27[) with 62 distributed uni- 
formly on the interval 0.005 < 62 < 0.165; these results were 
obtained by averaging (p n —po)/n (n = 120000) over an en- 
semble of 10 4 initial conditions (a;o,po) in the chaotic region, 
i.e., all having positive finite-time Lyapunov exponents. Solid 
line: The analytical result ()28|) . The inset shows the contin- 
uation of the main plot to smaller values of 62, distributed 
uniformly on the interval 0.00025 < b 2 < 0.00475. In this 
interval of very weak chaos, A is close to its maximal value 
1/4, see Eq. (p9)> . 
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FIG. 8: Circles: Numerical results for A versus a on the 
path 61 = 362 with 62 distributed uniformly on the interval 
0.005 < 62 < 0.1 (this path crosses both domains LR and R, 
see Fig. 3); these results were obtained as described in the 
caption of Fig. 7. Solid line: Analytical results from formulas 
and For small a, A w 0.12, close to the predicted 

limiting value of A = 0.125 from Eq. J30} (a = 1/3). 



available computational resources. An example is shown 
in Fig. 9. 

Our main result that the strongest Hamiltonian ratchet 
effect can arise in a limit of arbitrarily weak chaos has ap- 
parently no analog in ordinary ratchets if chaos is viewed 
as the deterministic counterpart of random noise. In fact, 
a sufficiently high level of noise is essential for the func- 
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0.11 0.14 



FIG. 9: Circles: Numerical values of A (obtained as described 
in the caption of Fig. 7) for K = 2.99 on the path b± = 0.25 
and 0.005 < &2 < 0.14; this path lies entirely in domain 7?, see 
Fig. 3 for the nearby value of K — 3. The Lyapunov exponent 
on the path varies in the interval 0.4787 < a < 0.8326. Solid 
line: Analytical results from formula (126[) for K = 3. 

tioning of ordinary ratchets or Brownian motors [l7l [l8j . 
Actually, it was recently shown that for a Levy ratchet 
the current decreases algebraically with the noise level 
(lij . in clear contrast with our results. 

The quantized version of our non-KAM system may 
be experimentally realized using, e.g., optical analogs as 
proposed in Ref. [IB] and is expected to exhibit in gen- 
eral a rich variety of quantum phenomena, including the 
quantum signatures of the weak-chaos ratchet accelera- 
tion. The study of these phenomena is planned to be the 
subject of future works. 
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APPENDIX A 

We derive here Eqs. (flO|) . First, continuity of the 
function ([§]) at x = b\ and x = 1 — b 2 implies that 

llbx=c — £>i, l 2 b 2 = 1 — b 2 - c. (31) 

Then, using Eqs. © and (|3"Tj) in the ratchet condition 
Jo f( x )dx = 0, we find that 

f f(x)dx = c-^c+^(l-c)-±=0, 

yielding the expression for c in Eqs. (|10j) . After inserting 
this expression in Eqs. (|3"TT) . we get the expressions for 
h and l 2 in Eqs. (TTOj) . 

APPENDIX B: THE REGION C FOR K = 3 



In this Appendix and in the next ones, we derive, for 
K = 3, exact results for the region C in Eq. (|2"3"|) . As 
mentioned in Sec. IV, several arguments and extensive 
numerical evidence indicate that C coincides with the 
chaotic region C for K — 3. We show here that one has 
the simple relation 

C = C = B U MB U M 2 B. (32) 

To show this, we first denote 

ijW = MB -BD MB, (33) 

B (2) = MB {1) -BC\ MB {1) . (34) 

We derive below the relation 

MB (2) C B. (35) 

Then, from the definition of C in Eq. (13"2"j) and from Eqs. 
([33 ]) - ((31) ]) it follows that 

MC C C'. (36) 

Eq. ([3"6"]) and the fact that M is area preserving imply 
that MC = C or M~ X G = C Thus, C = M S B U 
M S+1 BUM S+2 B for all integers s, which is possible only 
if C is equal to C in Eq. ([23]). Relation ([32]) is thus 
proven. 

To derive Eq. I|35p . we start by obtaining an explicit 
expression for B^ in Eq. (|33l) . For K — 3, the iterate 
of any initial condition (xq,Pq) under M satisfies p± = 
Xi — xq mod(l), Xi — xq + po + 3f(xo) mod(l). Clearly, 
when po varies in [0, 1) at fixed xq, x\ varies in the whole 
interval [0, 1). Then, taking (xq,Pq) in B and using Eqs. 
UHJl and ((33J, we get 

MB = {(x,p)\0 <x<l, x-h<p< x + 6 2 }mod(T 2 ), 

(37) 

B {1) = {(x,p)\h <x<l-b 2 , x-b!<p<x + b 2 }. 

(38) 

The region ([55]) is a strip (parallelogram) of slope 1, 
shown in Fig. 10(b) and corresponding to the strip B 
in Fig. 10(a). Next, we determine the region 

B (2) =MB {1) . (39) 

The second iterate (x 2 ,p 2 ) of (xo,Po) under M, with 
{xo,Po) € B and (xi,pi) G is given by 

p 2 = pi— 3(xi— c) mod(l) = — 2x\— x +3c mod(l), (40) 

x 2 = Xi + p 2 mod(l) = — 2xi + pi + 3c mod(l). (41) 
From Eq. (gDJ) and the first Eq. (jUJ), we find that 

p 2 = 2x 2 + x - 3c mod(l). (42) 
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which both the L and R accelerator-mode islands (AIs) 
exist (see Sees. IIIC and IV). For values of 61 and/or 




FIG. 11: (Color online) The region C, given by the union of 
the three regions shown in Fig. 10, with overlaps indicated 
by green (dark grey) color. Both AIs, L and R, exist in this 
case. See text for more details. 



FIG. 10: (Color online) The figure shows, for K = 3, 61 = 0.1, 
and &2 = 0.05: (a) the strip B, see Eq. {T3}; (b) the region 
B w , see Eq. or (|38 ]l : (c) the region B {2 \ see Eq. J39} 

or g3j. 

Eq. (|42"]) and the second Eq. (|4T|) imply that the region 
(13TJ1) is the following set of phase-space points: 

B {2) = {(x,p)} : 

\x = -2xi +pi +3c mod(l), (x 1 ,pi) 6 
I p = 2x + xq — 3c mod(l), —b 2 < xq < b\. 

The region (|43"]) is clearly a parallelogram of slope 2 folded 
into T 2 , as shown in Fig. 10(c). The region B^> in Eq. 
(f3"4"|) is given by Eq. ([4"3"]) with x restricted to the interval 
61 < x < 1 — 62- Then, using also Eq. (02), we get for 
(zcPo) G B and (x 2 ,J3 2 ) £ B^: 

MB (2) = {(x,p)} : 

lp = P-2 ~ 3(x 2 - c) mod(l) = x - x 2 mod(l), 
[ x = X2 +p mod(l) = xq mod(l). 

The result (JUJ) and Eq. dT3J) imply Eq. ([55]). 

APPENDIX C: AIs AND THE SHAPE OF C 

We study here the shape of the region C for K = 3 
in several cases. Let us write Eq. (|3"2"|) as C = C = 
_B U U i?' 2 ', i.e., the union of the three sets in Fig. 
10. This union is shown in Fig. 11, exhibiting a case in 



62 larger than those in Figs. 10 and 11, there may exist 
only one AI or no AIs (see Fig. 3). A case of C for 
which only the R AI exists is shown in Fig. 12 and is 
clearly different from that in Fig. 11. We show below 




0.2 f 0.4 0.6 0.8 1 
X 



FIG. 12: (Color online) The region C = B U U B (2) for 
K = 3, 61 = 0.3, and 62 = 0.03, with overlaps indicated by 
green (dark grey) color. Only the R AI exists in this case. 
See text for more details. 

that the existence of AIs and the shape of C in different 
cases depends on the location of the vertices (x^\p^) 
(j = 1,2,3,4) of the parallelogram (|4U)) . shown in Fig. 
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10(c), relative to the strip B. Because of Eq. (|39|) . one 
has (x^,p^) = M(xti),pW), where (#>'), pW) (j = 
1,2,3, 4) are the vertices of the parallelogram (|38|) in Fig. 
10(b). Clearly, 

= (6i,0), (* (2) ,P (2) ) = (b 1 ,b), (45) 
(x^\p^) = (l-fta.l), {x^\p^) = {l-b 2 ,l-b), 

where b = b\ + b 2 . To derive explicit expressions for 
(x^\p^) — M(x^\p^'), one has to properly determine 
the additive integers from the modulo operations in M 
so that (xti\pU>) will lie within the basic torus T 2 . We 
find that the values of x^> are 

x^=3c-2b 1 -l, x W=xV)=xW+b, 

X W = X W + 2b, (46) 

indeed satisfying < x^ < 1 in the relevant cases in 
which at least one AI exists. In fact, in these cases one 
has b < 0.5 [from Eq. ([211) with K = 3] and the latter 
inequality implies by simple algebra that the smallest 
value of iW in Eqs. (|46p. i.e., x^\ satisfies x^ > while 
the largest value (x^) satisfies x^ < 1. In addition, 
it is clear from Figs. 10-12 that the R AI exists only 
if x^ < 1 — 62 (vertex 3 is outside B); it is easy to 
show that the latter inequality is indeed equivalent to 
the existence condition 61 < F{b 2 ) for the R AI, derived 
in Sec. IIIC. Similarly, the L AI exists only if > 61 
(vertex 1 is outside B), which can be easily shown to be 
equivalent to the existence condition (|17[) . Thus, when 
both AIs exist, h < x& < 1 - 6 2 (j = 1,2,3,4). 

To determine the values of p^\ we first notice that 
the vertices (x^\p^') must touch the boundaries of the 
region ([3"T]) ; this is because the vertices (1431) in Fig. 10(b) 
obviously touch the boundaries of the strip B in Fig. 
10(a) and (x u \p^) = M{x^\p^). Then, in the case 
that both AIs exist, i.e., 61 < x^> < 1 — 62 (see above), 
(x^\p^') touch the boundaries p = x — b\ and p = x + b 2 
of the parallelogram (|38t [see Figs. 10(b), 10(c), and 11], 
so that 

p (l,2) = ^1,2) _ bu p (3,4) =x (3,4) +&2i (4?) 

Assume now that only the R AI exists, as in Fig. 12. 
Then, x' 1 ' < 61 (from above), i.e., vertex 1 (the point d 
in Figs. 11 and 12) lies within the left part of strip B, on 
the boundary of region (|37p given by p = x — &i mod(l); 
thus, for x^ < 61 (as in Fig. 12), pW in Eq. (g7]l must 
be replaced by x^ — b\ + 1 while p^ for j > 1 remains 
unchanged. Similarly, when only the L AI exists, vertex 
3 lies within the right part of strip B, on the boundary of 
region (|3"T)l given by p = x + b 2 mod(l); for a;^ 3 ) > 1 — b 2 , 
p(3) j n Eq_ (|47l) must be replaced by x^ +b 2 — 1. 

APPENDIX D: AREAS OF AIs 

Consider the L AI in Fig. 11. This is the triangle dfg 
on the torus T 2 , composed of two triangles, def and efg. 
The point d is vertex 1 in Fig. 10(c) and the segment de 



is part of the upper boundary of the region (|43|) . This 
boundary is a line of slope 2 passing through vertex 1: 

p-pW = 2{x-x^). (48) 

Then, since p e = 0, we get from Eqs. (|15 ]l -(|ig |l that 
x e = (3c — 61 — l)/2. Also, xj = bi and pf — 0. The 
point with x g = 61, lies on the line (|48)) with p^ 
replaced by p^ + 1. Thus, p g = 2 + 3&i — 3c. The area 
of the L AI is therefore 

S L = S def + S e fg = -{x e -xf)[p {1) + {l-p g )] 

= i(3c-3&!-l) 2 . (49) 

By symmetry arguments, the area of the R AI is obtained 
from Eq. (|49[) by inserting the expression for c from Eqs. 
(fTUl) and performing the exchange b\ ■n- b 2 . We get 

Sr = i(2-36 2 -3c) 2 . (50) 

APPENDIX E: AREA OF C 

The area of C can be calculated starting from the re- 
lation C = BU £?W U B (2) (see above), where £? (1) and 
B do not overlap. Then, because of Eqs. ([3"3"]) and (|3"9"|) . 
also B^> does not overlap with B^. However, it may 
overlap with B. The area of C is thus given by 

Sc = Sb + S'bu) + S B (2) — S BnB (2). (51) 

From Eq. (flU)) . Sb — b, where b = b\ + b 2 . The region 
£?W in Eq. (|3"8|) is a parallelogram with basis b (in the p 
direction) and height 1 — 6 (in the x direction), see also 
Figs. 11 and 12. Then, S Bm = 6(1 - 6). From Eq. 
and the fact that M is area preserving, it follows that 
S B (2) = S B (i). Finally, concerning the overlap B d B^ 2 \ 
we consider first the case that both AIs exist, see Fig. 
11. In this case, B n B^> consists of the green (dark 
grey) regions in Fig. 11. These are two parallelograms 
having heights b\ and 62 (in the x direction) and basis 
6, i.e., the width of region (|43|) in the p direction. Thus, 
S BnB (2) — b 2 . The area (|5ip is therefore 

5 c = 3(6 1 + 6 2 )- 3(6 1 + 6 2 ) 2 . (52) 

Consider now the case that only one AI exists, say the 
R AI as in Fig. 12. In this explained at the end 

of Appendix C, the point d, i.e., the vertex 1 of region 
B( 2 \ lies inside the left part of strip B, on the boundary 
of region (f3"T|) . This means that the black triangles def 
and efg in Fig. 12 are not included in the region B^ 
or B n B^> but they are actually part of the region (|37[) . 
Thus, to calculate 5' BnB (2) one must subtract from 6 2 
(the value of S BnB (2) in the previous case) the areas of 
def and efg. The area ([5Tj) is then obtained by adding 
Sdef + S e f g to the expression (|5"2"|) . By comparing Fig. 12 
with Fig. 11, it is clear that the areas Sdef and S e f g can 
be calculated precisely as in Appendix D and Sdef + S e f g 
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is given again by formula (l49l) . Therefore, the area (|52|) 
increases precisely by an amount equal to the area (|49[) 
of the missing L AI: 

S c = 3(6i + 62) - 3(6i + 6 2 ) 2 + (3c - 36i - if : /2. (53) 

Similarly, when only the L AI exists, one must add the 
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